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Abstract 

I propose a new version of the Hierarchical Reference Theory of 
liquids. Two formalisms, one in the grand canonical ensemble, the 
other in the framework of statistical field theory are given in parallel. 
In the latter the theory is an avatar of a new version of the non 
perturbative renormalization group (J. Phys. A : Math. Gen. 42, 
225004 (2009)). The flow of the Wilsonian action as well as that of the 
effective average action of Wetterich are derived and a simple relation 
between the two functionals is established. The standard Hierarchical 
Reference Theory for liquids (Adv. Phys. 44, 211 (1995)) is recovered 
for a sharp infra-red cut-off of the propagator 
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1 Introduction 

The first theory of liquids which incorporates the ideas of the renormal- 
ization group (RG) of Wilson [U [2] was proposed more than twenty-five years 
ago in a series of outstanding papers by Reatto, Parola and co-workers who 
developed the so-called Hierarchical Reference Theory (HRT) which gives an 
accurate description of the thermodynamics and the structure of a wide class 
of fluids close to and away from the critical region [3J, SI El El IZ] - The theory 
was developed and improved over years by these authors until very recently. 

During roughly the same period the Wilson approach to the (RG) has 
been the subject of a revival in both statistical physics and field theory. Two 
main formulations of the non perturbative renormalization group (NPRG) 
have been developed in parallel. In the first one, the RG transformations 
concern the Hamiltonian (or action) of the model [H [2j [HI [91 [10] while, in 
the second, the RG flow of the free energy is under scrutinity, or rather, to 
be more precise, the so-called "effective average action" [Til IT2"] . In both 
approaches, at scale "k" (in momentum space), a cut-off is introduced to 
separate the slow (q < k) and the fast (q > k) modes of the field. In the 
first approach the Wilsonian action is defined to be the effective action of 
the low energy modes yielding the same physics a long distances than the 
bare action. In this case the cut-off acts as an ultra-violet (UV) cut-off. In 
the formalism of Wetterich these ideas are implemented at the level of the 
Helmholtz free energy of the fast modes (in fact one is rather interested in its 
Legendre transform, the effective average action) and thus, in this approach, 
the cut-off acts conversely as an infra-rouge (IR) cut-off. The two points of 
view can be reconciliated formally [TOl US] . 

It turns out that many of new results obtained recently for the NPRG 
had been anticipated by the seminal works of Reatto 'et al.'. Since the two 
communities have few contacts and use a different language many things have 
been discovered or rediscovered independently by both sides in the ignorance 
of the results obtained by the other camp. Today field theorists would say 
that HRT is a sharp cut-off version of Wetterich effective average action 
approach to the RG, while specialists of liquids would claim that Wetterich 
theory is nothing but a field theoretical version of HRT with smooth cut-off 
and strange notations. 

A precise comparison between the two worlds however requires a field 
theoretical representation of simple classical liquids. Technically, a simple 
Hubbard-Stratonovich transform does the job [T31 H3] . This has been realized 
by many authors [16] H3 UHl EE] . The more achieved scalar field theoretical 
representation of liquids at equilibrium has been christened the KSSHE the- 
ory after the names of Kac [20] , Siegert [21] , Stratonovich [U] , Hubbard [15] 
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and Edwards [22] and is developed in references [TBI EE] inter alia. I show 
in the latter reference that the implementation of Wetterich effective average 
action method in the framework of KSSHE theory yields indeed HRT in the 
ultra-sharp cut-off limit. 

In general the field theory obtained via a Hubbard- Stratonovich transform 
yields a non-canonical theory in the sense that the coupling between the field 
<p and the source, here the chemical potential u, is non linear (see section 12.21 
below). An additional transformation is needed to obtain a standard field 
theory with a linear coupling. I adhered to this method in my previous works 
on liquids [181 US]- However it turns out that the NPRG ideas can easily be 
extended to the non-canonical theory itself and it has then the advantage to 
yield extremely simple functional relationships between the Wilsonian action 
and Wetterich's effective average action [13] . In the present paper I will show 
that applying this new field theoretical formalism to the KSSHE theory leads 
to a slightly different version of HRT with some advantages however. For 
instance the flow equations of Reatto 'et al' are reobtained by taking the 
limit of an ordinary sharp cut-off. Moreover the grand-canonical functionals 
corresponding to the Wilsonian and effective average actions of Wetterich 
are clearly identified which allows a deeper understanding of the NPRG for 
liquids. 

The paper is organized as follows. In next section [21 I circumscribe the 
type of liquids under study (the simplest!) and give their representations in 
the grand-canonical ensemble and in field theory. Then, in section [31 the 
coarse-graining operations of the various functionals of interest are detailed 
and their flow equations are derived in section [U in the last paragraph 14.31 
the usual HRT equation is obtained by taking the limit of a sharp cut-off. We 
summarize and conclude in section Two appendices explain some technical 
points. 

2 The model 

2.1 Grand-Canonical representation 

I consider a classical simple fluid made of identical hard spheres (HS) of 
diameter a with, in addition, isotropic pair interactions u(r^) (r^ = Xj\, 
Xi position of particle "«"). Since v(r) can be chosen arbitrarily in the core, 
we assume that v (r) has been regularized for < r < a in such a way that its 
Fourier transform v q is a well behaved function of q and that v(0) is a finite 
quantity. For convenience the domain Q C M d occupied by the molecules is 
chosen to be a d— dimensional cube of side L within periodic boundary (PB) 
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conditions; the volume of Q is noted V = L d . The fluid is at equilibrium in 
the grand canonical (GC) ensemble, /3 = is the inverse temperature 

(&b Boltzmann's constant), and p the chemical potential. In addition the 
particles are subject to an external potential ip{x) and we will denote by 
v(x) = (3(p — ip(x)) the dimensionless local chemical potential. I adhere to 
notations usually adopted in standard textbooks on liquids (see e.g. [23]) and 
thus denote by w(r) = —f3v(r) minus the dimensionless pair interaction. 
Moreover we restrict ourselves to the case of attractive interactions, i.e. 
such that w(q) > for all q. The model admits a thermodynamic limit if 
w(r) decays faster than \/r d+e as r — > oo, where d is the space dimensions 
and e > [21]. 

In a given GC configuration C = (N;x\ . . .xn) the microscopic density 
of particles at point x reads 

TV 

p(x\C) = Y,S {d) (x-x t ), (1) 

i=i 

and the grand canonical partition function (GCPF) H [u] which encodes all 
the physics of the model at equilibrium is defined as [23] 

E[u] = Tr[exp(-/W[C])] , 

oo 

Tr [■ • •] = E m / dl ■ ■ ■ dn e M-Pv m [C]] ■ ■ ■ . ( 2 ) 

where i = Xi and = c? d Xj. In equation (J2J /3V(hs) [C] denotes the HS 
contribution to the configurational integral. Its Boltzmann factor is either 
if configuration C involves overlaps of spheres, or 1 otherwise. Note that in 
fact S [v\ does not depend on the self energy — u>(0)/2. For a given volume V 
and a given inverse temperature /3, S [v] is a log-convex functional of the local 
chemical potential v(x) [25| 126]. I have employed here convenient matricial 
notations. 

V-p = / d d x z7(x)p(x|C) (3) 
Jn 

p-w-p = d d xd d yp(x\C)w(x,y)p(y\C). (4) 

2.2 Field theoretical representation 

It is possible to express H[V] as a functional integral making thus a link 
between the theory of liquids and statistical field theory [161 El EE] • With 
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the help of a Hubbard- Stratonovich transform [T31 US] one gets: 

SM = AC 1 / 2^ expf-^-i?-^ +W ( HS)^-w(0)/2 + ^n (5) 

where W(hs) = hiH(Hs), V 9 is a rea l random scalar field, R = w^ 1 is the 
inverse of w (in the sense of operators, i.e. w(l, 3) ■ R(3, 2) = 2)) and 

S(hs) I 27 ~~ w(0)/2 + ip] denotes the GCPF of bare hard spheres in the presence 
of a local chemical potential v{x) — w(0)/2 + (p(x). We optimistically suppose 
that the functional Sms) \y\ as well as its functional derivatives, i. e. the 
density-density correlation functions, are known exactly. This is nearly true 
since the HS fluid has been under scrutinity since so many years [23] . 
I have noted M w the normalization constant 

M w = j Vip exp ^-^v? • R ■ <fj ■ (6) 

The functional integrals which enter equations (EJ and fl6]) can be given a 
precise meaning in the case where the domain Q is a cube of side L with 
periodic boundary conditions, in this case the measure Dip reads [2] 

V<p ee T[-%fL (7a) 

d(pqd<p-q = 2 <i9R£> 9 for g 7^ . (7b) 

where A = (2tt/L) TL D [TL ring of integers) is the reciprocal cubic lattice. 
Note that the reality of <p implies that, for q ^ 0, £>_ 3 = where the star 
means complex conjugation. With this slick normalization of the functional 
measure one has exactly 

AC = exp J\n.w{qyj (8) 

in the limit of large systems (L — > oo). 

The field theoretical representation (jSJ) of the GCPF 5[i/] can be extended 
to arbitrary pair potentials. If w(l,2) contains a definite negative part cor- 
responding to a repulsive interaction then a second scalar field, purely imagi- 
nary must be introduced via an additional Hubbard- Stratonovich transform, 
see e.g. ref [18]. Handling these two fields makes the algebra a little bit 
more complicated but by no means intractable; we lazily skip this compli- 
cation in this paper and restrict ourselves to attractive w(r). Alternatively, 
one can gather the HS and repulsive pair potentials contributions in a single 
functional W^Ref) which will play the role devoted to W(hs)- 
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It is important to note that (jSJ) "looks" like the generator of Green's func- 
tion of a field theory with an action S [<p] = ^(p-w^-cp— W(hs) [v ~ w(0)/2 + <p\, 
with w(l, 2) playing the role of the bare propagator and the non local term 
Wms) that of the interaction. However purists will note that this is not a 
"canonical" field theory stricto sensu since the coupling between the field tp 
and the external source v is non-linear. As detailed at length in refs. [121 [32] 
this point can be circumvented by the change of variables x — <fi + v — w(0)/2. 
By doing so one recovers a "canonical" field theory involving a linear cou- 
pling between \ and v. Moreover one can establish a rigorous mapping 
between these two field theories. Then all the knowledge amassed in sta- 
tistical field theory can be injected in the theory of liquids, notably well 
established perturbative technics and/or more recent non-perturbative ap- 
proaches. Recently however we proposed rather to consider the action S[<p] 
as a well-educated one and work directly with it rather than to try to schol- 
arly adhere to a canonical action via the mapping </?•<->■ x- As discussed 
in reference [13], the adoption of this "non- canonical" point of view leads 
to a remarkably simple reconciliation between the points of view of Wilson- 
Polchinski on the one hand and that of Wetterich in the other hand. Applied 
to liquids, as will be done here, it yields a new version of the smooth HRT 
theory with some advantages on my previous attempts [19J. 

3 Coarse graining 

The coarse graining procedure involved to build families of actions with 
the same physics at long distances is the core of the RG theory. In this section 
we introduce this procedure step-by -step following recent developments in 
field theory [El [13] 

3.1 Ultra-violet regularization of the pair potential 

I noted in section I27T1 that the pair potential (propagator) w(r) can be 
regularized at short distances < r < a, i. e. in the ultra-violet (UV) regime, 
without changing the GCPF of the system. I introduce a potential Wq(t) 
defined in Fourier space by 

w£(q) = C(q,A)w(q) , (9) 

where C (q, A) is an UV cut-off which is equal to 1 for q < A and to for 
q > A. Typically C(q,A) = 1 — Q e (q,A) where e (g,A) is a function which 
looks like the step function 0(g — A). It could be precisely this function 
and we would then speak of a sharp cut-off; or it could be a smooth version 
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with Q e (q, A) varying smoothly from to 1 in a small interval (A — e, A + e) 
about the UV cut-off; in this case we would speak of a smooth UV cut-off. 
It is sometimes convenient to choose C(q,A) = C(x = q/A) but by no means 
compulsory. 

In any case if Act 3> 1 then w(r) and w$ (r) coincide outside the core 
and moreover we have Wq(0) = w(0) and \wq (r = 0)| < oo. Although this 
UV regularization does not change the GCPF H [u] we shall henceforth note 
this functional Sq [is] to emphasize that it depends functionally on Wq. For 
further convenience, we also need a special notation for the inverse of the 
regularized propagator which will be christened Rq = [wq^ 1 . We can 
thus rewrite (j2J) and (jSJ) as 



Of course 



Tr 



exp ( lp. w £.p + p. [ v -v4(0)/2] 



(10a) 



j^Jvip exp \^yp -B%-<p+ W { ns) [v - w^)/2 + ^0b) 



Sq [u] = S [z/] as soon as Act ^> 1 which will be assumed. 



3.2 Two lemma 

I prefer, for the convenience of the reader, to give here two lemma nec- 
essary for subsequent developments rather than to postpone them in an ap- 
pendix. Let A some definite positive operator and J-'lf] an arbitrary func- 
tional of real scalar field <p, then 

If 1 
Jf^J V ¥ ex P(-2^ • A_1 • <P) + Vol = exp(A))^M , (11) 

where the operator Dq reads 

D Q ... = \f d d xd d yA(x,y)—-f——-. (12) 

An interesting application of lemma (fTTj) is to reconsider equation ( llObj) as 

Efa + w£(0)M = e D °Z m [v] , (13) 

therefore the operator e D ° constructs the full GCPF of the system from the 
GCPF of the reference HS system. 

The second lemma is sometimes referred to as Bogolioubov theorem : 
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Let Aj, 1 < i < n, be n definite positive operators and J 7 ^] an arbitrary 
functional of real scalar field <p, then 



n{^exp(-^.A-.^)U ( f: 

i=l ^ ~ ^ i=l 



Pi 



V( P , 1 A-l % 

_exp(--<p.A -^)x 



x ^"(v) 



(14) 



where A = Ya=i ^i- 

The two lemma ffTTl) and (|14|) are easy consequences of Wick's theorem. 



3.3 Blocking the action 

I now apply the exact RG approach of Tim Morris [TUl [T5] to our non- 
canonical field theory. As a consequence of Bogolioubov theorem (j!4p the 
GCPF Hq [u] can be rewritten in terms of two propagators and two fields as 



J V V< exp (- l - V< ■ R k ■ V< ^j ~£ [u - w k (0)/2, V< ] 



4 



Mr 



+ w< 



(HS) 



Vify exp ( —£P> - R k-P> + 



where < k < A is the running scale of the RG and where 



ip = if < + ip > and Wq = w k + Wq . 



(15a) 



(15b) 



(16) 



In (|T5l) - (|T6j) I have separated the field <p into "rapid" (</>>) and "slow" modes 
(</?<)■ The low-energy modes are associated to the propagator Wq (with 
inverse Rq) which is cut off from above by k, while the high-energy modes are 
associated to the propagator w£ (with inverse R^) which is cut off from below 
by k and from above by A. I demand that w^(q) = wo(q)(C(q, A) — C(q, k)) 
should be positive which will be assumed henceforth. In the popular case 
where C(q,A) = C(q/A) it is sufficient for the cut-off function C(x) to be 
a monotonous decreasing function of its argument. Some comments are in 
order. 

(i) In order to establish equations ffl~5|) I used Wq (0) = w£(0) + Wq(0) as 
implied by ffT6l) . 
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(ii) Note that the "slow" and "rapid" modes, respec. (p < and <£>> are or- 
dinary scalar fields, in particular they have a full spectrum of Fourier 
components <p<(q) and <p>(q) (even the condition < q < A is not 
compulsory if the bare propagator is UV regularized). The cut-off at 
scale "k" acts only on the propagators. 

(iii) The functional S^fv, tp\ noted in this way by Morris [TU] is the crux 
of the whole matter since it allows, as I will show, to make explicit 
the link between the Wilsonian action and the effective average action. 
However here this link proves trivial since, as apparent in formula (115bl) . 
E£[v, ip], depends functionally only on the sole variable v + <p. 

Let me first set y?< = in (|15bp . It yields 



~£ [u, V< = 0] 4 s£ [u] (4 exp {W k A M) ) (17a) 
= -±-Jv<p exp {-lip -R^.ip + W m) [u - <(0)/2 + ^ 

(17b) 

From the point of view of field theory this shows shows that W k H = 
In [v] is the Helmholtz free energy of the rapid modes <p> in the presence of 
the source u(x); W k is thus the generating functional of connected correlation 
functions with UV regularization (at A) and an infra-red (IR) cut-off (i. e. 
at k). This functional is related by a Legendre transformation to the effective 
average action of Wetterich, as it willl be discussed in section SJ 

From the point of view of the theory of liquids clearly, (see e. g. equa- 
tion (jSJ)) [u] is precisely the GCPF of a system of hard spheres with 
additional pairwise potentials w£ [r] , i. e. , the k-system, to adopt the 
terminology of Reatto 'et al.' Therefore one also has : 

i— 'A 



S£ \iA = Tr 



exp[-p-w£-p + p-[u-w£(0)/2] 



Note that equation ffl8|) is valid provided the propagator w^(x, y) (pair po- 
tential) is definite positive (attractive) so that a Hubbard-Stratonovich trans- 
form is licit; that is why I imposed earlier that the cut-off function x — > C(x) 
should be a decreasing function. Comparing equations (115bj) and (117bj) one 
thus has synthetically 

E# [i/, d = Sf [i/ + d . (19) 
I now introduce the Wilsonian action 



(20) 
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which allows me to rewrite ( I15al) as 

"0 [v + 4(0)/2] = j^- Jv V< exp f-^< • R k ■ - S£[v + </?<] .) 

(21) 

Indeed, in equation ( 121]) S 1 ^ [</?<] play the role of the effective action of the 
slow modes at scale k [TJ [2j [10]. Here k plays the role of an UV cut-off. 
Note that the functional identity Sfc = — is not true for a canonical field 
theory [TO] . 

We end this section by applying lemma ( TIT]) to the 2 equations ( TT5]) which 
can thus be rewritten as 

3fri/ + «;*(0)/2] = e^SfM (22a) 
^[^ + <(0)/2] = e^5 (HS) H (22b) 

from which it follows that 

S A [z, + ^(0)/2] = e B oe^S { Hs)M. (23) 
Comparing with equation ( TT3"]) yields 

e^ A = e D o e ^ . (24) 
This is the semi-group law of the RG. 



4 RG Flow equations 

4.1 RG flow of the grand potential W£ 

From the definitions ( [16]) when A; — > A then w£(r) — >■ and the k-system 
at = A is a fluid of hard spheres. When k decreases from A to k = 
more and more Fourier components are included in the potential w£(r) and 
finally precisely at k = the k-system coincides with the full model since 
tffc(r) tends to Wq{t), i. e. essentially w(r) since both potentials differ only 
in the hard cores with no effect on the physics. In this section I establish 
the equations which govern the flow of the grand potential W- introduced in 
Sec [3] and its Legendre transform. The flow equations of the Wilsonian action 
<S^ follows trivially. Note that the formal solution of the not yet established 
flow equation for Wfc is already known, it is given by equation ( ]22bj) . 

The flow equation for can be obtained in the framework of statistical 
field theory and, apart some tricks due to the self-energies, along essentially 
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the same lines as in reference [IB]. I postpone this derivation to appendix 
and I give here only the proof in the framework of liquid theory. An easy 
task. The starting point is equation (fl8l) which I differentiate with respect 
to k at fixed u(x). It gives readily 

dkW£[v]\ u = l - I d d xd d yd k wt(x,y) {Gt(x,y) - pi(x)5^(x -y)} (25) 

where p k {x) =< p(x\C) >gc denotes the mean density of the k-system in 
the GC ensemble and G k (x,y) =< p(x\C) p(y\C) >gc its pair correlation 
function at chemical potential v(x). Since one has p k (x) = 5W£ / 5v(x) = 
Wj»\x) and also G k T)A (x,y) = G k T) ' A (x,y)-p£(x)p£(y) = 5^W k A /8u(x)8u(y) 4 
Wj, A (x,y) one can rewrite equation (125]) under the closed form 

d k W£[v + w£(0)/2]\ v = \fd d xd d yd k w K {x,y){wi 2) \x,y) 

+ \x)wi l) \y)} . (26) 

This partial differential equation (PDE) must be supplemented by the 
initial condition W k = Wpas)- This equation is closed only superficially 
since the pair correlation function W^ A (x) and A (x,y) depends func- 
tionally on the chemical potential v{x). Differentiating functionally succes- 
sively both members of equation ( 126"]) with respect to the field v(x) one 
obtains a hierarchy, or tower, of equations for the A (x\, . . . ,x n ) = 
8 n W k K /8v{x{) . . .8v{x v ), i. e. the Green's -or connected correlation- func- 
tions. 

4.2 RG flow of the Wilsonian action 

The flow equation ( 12"B"|) has the structure of Wilson-Polchinski's equation 
for the Wilsonian action; indeed with S k = — W k it takes the usual form 

(BUSHED! 

d k S£& + w£(0)/2]\ (p = \ j d d xd d yd k wUx,y){S { k 2)A (x,y) 

- S^ K {x)Sf K {y)). (27) 

where A (xi, . . . , x n ) = 8 n S£/8(p(xx) . . . 8ip(x n ). This equation, which 
must be supplemented with the initial condition = — W(hs)> does not 
appear to have yet been considered in the theory of liquids contrary to the 
abundant literature devoted to it in statistical field theory (see e. g. [9] and 
references quoted herein). 
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4.3 RG flow of the effective average Kohn-Sham free 
energy f3A% 

Since the seminal works of Reatto 'et al.' [31 IU El E] it is usual to consider 
rather the flow for the Kohn-Sham free energy A k [p] of the k-system which 
is defined as the Legendre transform of W£ [u] . One has the usual couple of 
relations 



sup <^ v ■ p - pA k [p] \ 



PA$\p] = sup [p-P-W£[v]} Vp 



(28a) 
(28b) 



n-rA 



Obviously the flow equation for j3A k [p] should be deduced from that of 
W k [u]. Let me consider for instance equation (128aj) . I have for all v 



,-rA 



W k A [u]=p*-v- &A k [p 



(29) 



where p*(x) is, if it does exists, the unique solution of the implicit stationnary 
equation 



waI [ P ] 



5p(x) 



(30) 



Differentiating ( 129]) at fixed v gives 



,-rA 



d k W£[v\l = d k p*-v- 



50A k [p 



5p* 



d k p* - d k /3At [p* 



(31) 



which further simplifies thanks to the stationnary condition (I30p with the 
final result : 



dkW k A [u] | y . 



(32) 



Similarly, starting instead from equation f l28b|) one finds that 



-rA 



dkA k [p] 



Ar, .*l 



- d k W k A [u 



(33) 



where v* is the unique chemical potential which is, if it exists, solution of 
equation (128b|) . I shall drop the subscript "V in the sequel. 
At this point I introduce the direct correlation functions 



7?WA/ 

C k [xi, 



,-rA 



x n ) = S n A k/Sp(xx) . . . 5p(x n ) 



(34) 
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which are the analogous to the vertex functions in field theory. I have simpli- 
fied the notations and it must be stressed that, in addition to their arguments 
these function depend functionally on the profile p(x). The C k A {x\, . . . , x n ) 
and the W^ A ( linked by generalized Ornstein-Zernike (OZ) 

equations and in particular one has the usual OZ equation : W^ A (x,y) = 
— [C k A ]" 1 (x, y) [6J. This last property allows to rewrite equation (1321) in a 
closed form involving A k , its functional derivatives and the propagator. One 
finds : 



d k A A k [p] =-\f d d xd d yd k w k {x,y)[C k 2)A }-\x,y) 
p 1 Ju 



+ - / d d xd d yd k w k (x,y)p(x)p(y) -\ [ d d x d k w%(x, x)p{x) . 
1 Jn 1 Jn 

(35) 

Note that, in order to obtain the result I made use of the trick dkiv^x, y) = 
—dkW A (x,y). However, the resulting equation f )35|) is quite awkward and it 
is convenient, by adapting the ideas of Reatto 'et a!.' [21 IH El E] to our case, 
to introduce a modified Kohn-Sham free energy : 

(3At [p] = 0A A k [p] -\p-wlp+ l -p- w k (0) . (36) 

This new functional was introduced by Wetterich 'et al.' [TTJ in statistical 
field theory under the name of "the effective average action" independently 
from Reatto 'et al.'; it differs from (3A k [p] by a simple quadratic form and 
satisfies obviously the following flow equation 

d k Ai[p] = - 1 - f d d xd d yd k w k (x,y) {c { k 2)A -w k o y\x,m a ) 

An) As \ $ n A k lp] . ,s 

CI' [xi,...,x n ) = x 7 \ 37b 

I already obtained this equation in my first version of the RG group for 
liquids but with the help of a canonical field theory [TU]. Note that the direct 

correlation functions C^ A (x\, . . . , x n ) differ from the C k A {x\, . . . ,x n ) only 
for n > 3; in particular one has 

C ( k 2)A (x,y)=C { k )A (x,y)+w k (x,y). (38) 

I stress that although f3A k [p] is a convex functional of the profile p(x) this is 
not the case in general for (3A k [p] except of course at k = since (3A k=0 [p\ = 
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j3A A =0 [p\ as apparent in formula fl3"ET) . Indeed although C^f A < in the sense 
of operators, equation (155]) shows that could become positive due to 

the addition of the positive operator w^. 

Equation f)37ap must be supplemented with an initial condition at k = A. 
From W£ = W(hs) as follows from equation (I22bj) and D 1 ^ = one concludes 



rA 



that (3A A = /3A(h_s) an d thus from 

MJW = /?Ahs)[p] - • ^ A • P + • <(0) • (39) 

In KSSHE theory the initial conditions are identical but are obtained in a 
more complicated way which requires the divergence of the regulator at k = A 

It turns out that the initial Kohn-Sham free energy of the k-system, i. e. 
at k = A, coincides with the mean-field Kohn-Sham free energy at any scale 

"k", i. e. one has /3^ A(MF) [p] = Mjy[p] for a11 k - Moreover this MF free 
energy constitutes a rigorous upper bound to the exact free energy and thus 
/3^[p] > f3A A [p}. This result already obtained in references [TBJ [19] for the 
KSSHE theory is proved for this new version of the theory in appendix 1X1 
where a precise definition of the MF approximation is provided. 

For a homogeneous system (3A^ [p] = Vfj}(p) where (3 times the free 
energy per unit volume f^(p) is a function (not a functional) of p and the 
flow equation for /■ (p) is easily deduced from equation f)37al) and reads 

Cfc/A - -n / ~( 2 )A, x ^ ' ^ 4U ) 

where J = J d d q/ {2ir) d and C^ A (q) is the Fourier transform of A (c, = 

C^ A (r = \\x — y\\) for a translationally invariant fluid. Note that the right 
hand side is non singular since the denominator is negative definite. It is a 
general property that the flow of f£ has neither UV or IR singularities, in 
particular IR singularities (near a critical point) are smoothened by k and 
they build up progressively as the scale-k is lowered [TT| [12]. . 

To make contact with the work of Reatto 'et al.' let me consider now the 
case of a sharp cut-off, i. e. C(q, k) — 1 — Q(q — k). Then equation (l4"Uj) 
becomes [T5] . 

** - ^> {1 -#^ (41) 

where Sd = 2ir d ^ 2 /T(d/2). This is the usual HRT equation [6]. In order 
to obtain fj4Tj) from (j4l?l) one replaces the sharp Q(q — k) by a smooth one, 
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Q e (q,k), denotes S e (q,k) = — <9fcO e (g, k) the smooth 5-function and finally 
make use of Morris lemma which states that, for e — > 



provided that the function f(Q e (q, k), k) is continuous at k = q in the limit 
e — > 0, which is the case here. In the canonical field KSSHE theory for liquids 
one needs to consider rather an ultra-sharp cut-off regulator [13J. 

Whatever the considered version of HRT, smooth or sharp cut-off, i. 
e. respec. equations (j4T)]) and (j4"T|) . both equations are of a formidable 
complexity despite their apparent simplicity Indeed the kernels of these PDE 
involve the two-body C^ A [p] which depends functionally on the density. 
Clearly one can deduce from these equations an infinite tower of equations 
for the Cfc A [p] (1, . . . , n) by differentiating them functionnaly with respect 
to the density profile p(x). A discussion of this point is however out the scope 
of the present paper, it has been developed in detail in a general context in 
reference [15] . 

5 Conclusion 

In this paper I have developed a new version of the "smooth" hierarchical 
reference theory for liquids. It can be seen as an effort to conciliate the point 
of view of the theory of liquids and that of statistical field theory. Technically 
it is an application to a peculiar field theory (KSSHE), aimed at representing 
liquids, of the more general theory developed in reference [13]. The equiv- 
alence pair potentials propagators is the crux of the whole matter. In 
addition, the interplay between the points of view of Wilson and Wetterich 
in their application to the theory of liquids yields a better understanding of 
the work of Reatto and his collaborators. 

From a prosaic point of view the differences with the first version of 
this work [19] are small but with interesting consequences. In practice, if one 
leaves aside philosophical considerations on Wilsonian and Wetterich actions, 
the main changes concern essentially the expression of the pair potential w£ 
of the k-system. In the first version, following Wetterich, I add a mass term 
<p ■ K k ■ ip/2 to the KSSHE action. With K k (q) ~ Zk 2 (l - G e (q,k)), Z 
large. This results in an ugly pair potential for the k-system, of the form 
^k(l) = ™o (<?)/(! + ( ( l)'R-k(q)) ■ In order to recover Reatto's equation (jlTp 
one then needs to consider an ultra-sharp cut-off limit, i. e. : e — > and 
Z — > oo. I recall that, here, only the simple sharp cut-off limit e — > 0, 
[Z = 1) is required. 




(42) 
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In the present formalism well-chosen cut-off functions such as C(x) = 
exp(— x 2 /2) yields analytical w£(r) in many cases, for instance w(r) = exp(— ar)/r 
(Yukawa), w(r) = l/r d+n , etc. Indeed one then obtains for w£(r) = ifEwaidtj) 
nothing but the Ewald potential in direct space, the guy of numerical simula- 
tions. The lacking contribution in Fourier space corresponds to interactions 
between "blocks" of size 1/k. This potential iUEwaid( r ) is monotonous in gen- 
eral. In the old version the potential w^{r) is not analytical and oscillates in 
direct space. 

A last difference can be emphasized : in the old version the initial condi- 
tion for the flow equation of (3A^ requires that the cut-off function 71 a = oo, 
which is difficult to implement rigorously in practice and could lead to unex- 
pected errors in numerical applications. In the present version, as discussed 
in section 14.31 there is no such problem and the initial condition is easily 
implemented. 

Despite the beauty and success of Reatto's HRT some drawbacks of the 
theory appeal a smooth cut-off version of it, such as the one exposed here. 
A sharp cut-off regulator induces singularities in the pair potential wfc(r) 
which decays slowly as ~ cos(kr)/r for r — > oo. This circumstance makes it 
difficult to solve exactly the integral equations involved in the most widely 
used closures of the Hierarchy; one has to resort to approximations which, 
as discussed by Reiner [27] are sometimes not easy to controll. A second 
drawback of the sharp cut-off version is the value of some critical exponents, 
notably that of the specific heat, a ~ —0.05, which is negative, while positive 
as it should in the smooth cut-off version [28J . 

It seems to me that solving exactly the smooth cut-off version of HRT is 
possible; that is coupling the numerical solution of a PDE (flow equation) 
and integral equations (closure of the hierarchy). Some of the drawbacks of 
the Reatto's 'et al' version of HRT could then be cured. 
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A Alternative derivation of the RG flow of 

I derive the flow of W k in the framework of field theory (see also reference [13 
I will start from equation (122b|) 



L = R 

Et[u k = u + wt(0)/2] = e^H (H s)M, (43) 

and take the partial derivatives of the both sides "L" and "i?" of EQ. (143]) 
with respect to the scale "k". One thus has for the l.h.s. 

5Z£ 



d k L\ u = d k E%[is k ] + I d d x 6 -f^d k u k 



H^N x {d k W£[v k \ + i / d d xW<; 1)A [v k ](x)dw£(0)/2} , ( II) 



2 

and, for the r.h.s. 

d k R\ v = d k D£e D £z (HS) [v 



n 



' f d d xd d yd k w A (x,y) 



2 Jn ' 5is k (x)5v k (y) 

1 
2 



x I I d d xd d yd k w£(x,y)G£iv k }(x,y) , ( lo) 



n 



where I used the expression : 

Dt--- = \l d d xd d yw£(x,y) f ^ . (46) 

Equating equations (jUJ) and and noting that the equality, obtained for 
arbitrary u k , is then valid in the change v k — > v arbitrary, one obtains 
finally 

dkWM\ v = \ I d d xd d yd k w£(x,y) {G£(x,y) - pt(x)6( d) (x -y)} , 

(47) 

which indeed coincides with equation fl25]h 



B Mean field theory at scale "k" 

The Mean-Field approximation of the GCPF of the k-system will be de- 
fined as 



r-nA (MF) 



M = exp {-~y>* • J# • <p* + W m) H + ¥»*]} 



(48 
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where v k = v — w£(0)/2 and <p* is the location of the saddle point inte- 
grandl f!17bj) . <p* satisfies the implicit relation 

R*-<p* = -g^w m [u k + <p*] 

= PQ3S)Wk + <P*\(?). (49) 

For a given chemical potential v{x) the MF profile of the k-system is given 

by p^ M (x) = 51nH^ MF ^[i/]/5i/(x) = P(as)[ u k + L P*]{ X ) where I made use of 
the stationary condition ( l4"9"j) . Therefore the "true" Kohn-Sham free energy 
of the k-system is given by 

/3^ (MF) = -lnS^ MF) + z,.pf MF) . (50) 
A short calculation will show that [191 H3] 

/?X (MF) [p] = MchsjM -\p-<-p + \p- <(0) • (51) 

where /3Ams) [p] is the free energy functional of the HS fluid at same density. 
Therefore one finds for the average effective Kohn-Sham free energy 

/3^ (MF) [p] = /3Ahs)H - \ P ■ < ■ P + \ P ■ w£(0) . (52) 

Note that f3A^ ^ M [p] is independent of scale "k" and thus equal to its initial 
value PA A A iMF) [p]. 

I prove now that 0A^ MF [p] is a rigorous upper bound to the Kohn- 



Sham free energy of the k-system [26| [18] . We consider the Legendre-Fenchel 
relation for the reference system of hard spheres : 

W(hs) [v] = sup {v ■ p - PAqjs) [p] } ( Vi/ ) (53) 
p 

which implies the Young inequalities 

( Vi/ , Vp ) W m [u] + 0A m [p]>p-v- (54) 
Injecting this inequality in the espression (117bj) of [u] one gets ( W, Vp) 

jj^ exp(--^ • R% ■ cp + p ■ cp) 
> exp(-M(HS) [p] + P ■ + ^P ■ w£ ■ p) , (55) 
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where I made use of Wick's theorem. Rearranging the terms one thus has 
(Vz/.Vp) p-v-W£[v\<fa K k m) [ P \, (56) 

(Vp) 0A A k (MF) [p] > sup{p ■ v - W£ [u] } . (57) 



Hence, from the very definition of the Legendre transform f!28bj) 



(Vp) /^ (MF) [p]>/3^[p] ( 58 ) 

Turning now our attention to the effective average Kohn-Sham free energy I 
find the exact upper bound : 

(Vp) PAt< m \p]=PAi\p]>/3Ai\p]. (59) 
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